Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN M?eX style file vl.4) 



Effects of dynamical evolution on the distribution of 
substructures 



Jorge Penarrubia^ , Andrew J. Benson^ 

^ Max Planck Institute fiir Astronomie, Kdnigstuhl 17, D- 69 117, Germany 
^Department of Physics, University of Oxford, Keble Road, Oxford, OXl 3RH, U.K. 



o 
o 

(N 
o 

Q 

in 



> 
o 

m 
(N 

O 

o 



2 February 2008 



ABSTRACT 

We develop a semi-analytical model that determines the evolution of the mass, position 
and internal structure of dark matter substructures orbiting in dark matter haloes. We 
apply this model to the case of the Milky Way. We focus in particular on the effects 
of mass loss, dynamical friction and substructure-substructure interactions, the last 
of which has previously been ignored in analytic models of substructure evolution. 
Our semi-analytical treatment reproduces both the spatial distribution of substruc- 
tures and their mass function as obtained from the most recent N-body cosmological 
calculations of Gao et al. (2004). We find that, if mass loss is taken into account, the 
present distribution of substructures is practically insensitive to dynamical friction 
and scatterings from other substructures. Implementing these phenomena leads to a 
slight increase (~ 5%) in the number of substructures at r < 0.35rvir, whereas their 
effects on the mass function are negligible. We find that mass loss processes lead to 
the disruption of substructures before dynamical friction and gravitational scattering 
can significantly alter their orbits. 

Our results suggest that the present substructure distribution at r > G.35rvir reflects 
the orbital properties at infall and is, therefore, purely determined by the dark matter 
environment around the host halo and has not been altered by dynamical evolution. 

Key words: method: semi-analytical calculations - galaxies: kinematics and dynamics 
- cosmology: dark matter 



1 INTRODUCTION 

The cold dark matter paradigm for the formation of struc- 
ture in the Universe has met with considerable success in 
describing both the large scale distribution of matter in our 
Universe (Bacon, Refregier & Ellis 2000, Percival et al. 2001, 
Spergel et al. 2003) and the structure of galaxies and galaxy 
clusters (see for example Klypin, Zhao & Somerville, 2002, 
Kneib et al. 2003; although potential discrepancies remain, 
see for example Weldrake, de Blok & Walter 2003, Sand et al. 
2004). In recent years, high resolution N-body simulations 
have demonstrated that a record of the hierarchical merg- 
ing process through which dark matter haloes form in such 
cosmologies is preserved in the form of substructure within 
haloes (Moore et al. 1999, Klypin et al. 1999, Springel et 
al. 2001). When a small halo merges into a large, "host", 
system it is not immediately destroyed, but instead begins 
to orbit within the host, gradually losing mass due to the 
actions of tidal forces. 

These substructures are presumably the locations of 
galaxies in clusters and of dwarf galaxies in the Local Group 



(Kauffmann, White & Guiderdoni 1993, Bullock, Kravtsov 
& Weinberg 2000, Somerville 2002, Benson et al. 2002b). 
The properties and evolution of substructures have been ex- 
tensively studied in recent years. The distribution of halo 
masses and the distribution of substructures throughout 
the host are now well studied (though that is not to say 
agreed upon) through both N-body and semi-analytical 
means (Ghigna et al. 2000, de Lucia et al. 2004, Gao et 
al. 2004, Taylor & Babul 2004). 



Although substantial computational efforts have been 
made to analyse the present distribution of substructures 
in dark matter haloes, the subject still lacks detailed the- 
oretical studies that explore the main processes that drive 
the evolution of those systems. In particular, one aspect of 
substructure evolution that has not been examined in de- 
tail is the interactions between substructures. Although such 
interactions are self-consistently included in N-body simu- 
lations their importance has not been assessed, while pre- 
vious semi-analytical calculations have ignored interactions 
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between substructures*. These encounters are an additional 
source of heating for substructures, thereby altering their 
rate of mass loss, and also change the distribution of sub- 
structure orbits (and therefore their spatial distribution). 

In this work we examine whether dynamical evolution 
occurring within the host halo alters the present distribu- 
tion of dark matter substructures. We use a semi-analytic 
model that implements the main dynamical processes acting 
on substructures: dynamical friction, gravitational scatter- 
ing, gravitational heating and mass loss. Semi-analytic codes 
have several advantages over N-body calculations: (i) they 
allow us to isolate the effects of each phenomena, (ii) they 
do not suffer from resolution limitations and (iii) they re- 
quire fewer computational resources. However, these codes 
are based on analytic approximations to complex dynam- 
ical processes. This inevitably limits the accuracy of such 
codes. However, such models have been shown to provide a 
good description for substructures in cosmological dark mat- 
ter haloes (Taylor & Babul 2004, Benson et al. 2002a). We 
will quantify how the neglect of substructure-substructure 
interactions impacts upon the results of such models. 

Ultimately, the question that we are attempting to an- 
swer in this work is the following: Is the present distribu- 
tion of substructures driven by physical processes acting in 
the host halo (in the sense that the orbits of those systems 
are strongly altered by dynamical processes), or, alterna- 
tively, does their distribution reflect their orbital properties 
at the time of accretion? The answer to that question is, for 
instance, of great importance for observational surveys of 
satellite galaxies, which use these bodies as a direct indica- 
tor of the properties of the host haloes and the dark matter 
cnvirornncnt (e.g. Prada ct al. 2003, Sales & Lambas 2004). 
In that sense, knowing whether dynamical effects are impor- 
tant for the substructure distribution could help understand 
whether such effects can significantly alter the properties of 
observed satellite galaxies. 

The remainder of this paper is arranged as follows. In 
§2 we describe how we generate host haloes and their popu- 
lation of infalling substructures. In §3 we describe the semi- 
analytic code used to evolve the substructures. In §4 we 
describe the calculations carried out in this work. We first 
show an illustrative study of the substructure-substructure 
encounters that these systems suffer during their evolution. 
We then analyse the evolution of the density profile and the 
orbit of substructures and examine how dynamical friction, 
gravitational scattering and mass loss affect the present spa- 
tial and mass distribution of substructures. We also examine 
the evolution of the substructure mass function. Lastly, in 
§5 we summarise our results. 



* Springcl ct al. (2001) included mergers between sub-haloes in 
their combined N-body-scmi-analytic study of a z = galaxy 
cluster, but did not assess in detail the importance of these merg- 
ers. Somerville & Primack (1999) also allowed for the possibility 
of mergers between substructures in their semi-analytic model of 
galaxy formation. Note that neither of these studies addressed 
the questions we are interested in here, namely the effects of 
substructure-substructure encounters as opposed to mergers. 



2 INITIAL SAMPLES OF DARK MATTER 
SUBSTRUCTURES 

We generate initial conditions for our semi-analytical cal- 
culations by constructing merger trees for a sample of dark 
matter haloes identified at z = 0. Specifically, we begin with 
dark matter haloes of mass 2x 10^^h~^ Mq, similar to that of 
the Milky Way (Klypin, Zhao & Somerville 2002). We then 
apply the methods of Cole et al. (2000) to construct the 
merging history of this halo back to high redshift. For this 
calculation we assume a ACDM cosmology with parameters 
Qo = 0.3, Ao = 0.7, 0-8 = 0.9, _ffo/100km s-^ Mpc^^) = 
0.7 and F = 0.21, together with an inflationary initial power 
spectrum (i.e. P{k) oc A;" with n = 1). At each point in time, 
we identify the most massive progenitor of the z = halo 
and classify this as the "host" halo at that redshift. We then 
identify haloes which are about to merge with the host at 
each time and label these as new satellites. 

At each point in time, we record the virial mass, virial 
velocity and concentration of the host halo (which is as- 
sumed to have an NFW density profile). We also record the 
same information for each satellite along with the epoch of 
each merging event. We select an initial orbit for each satel- 
lite according to the method of Benson (2004) — specifically, 
the initial position is chosen at random on the surface of 
a sphere with radius equal to the current virial radius of 
the host halo, while the initial velocity is chosen from the 
distribution found by Benson (2004) for 2; = in ACDM 
cosmologies^ . 

In this way wo generate both the formation history of 
the host halo and the properties of the infalling population 
of satellites. These are used as inputs to our model for the 
evolution of substructures. 

We generate a sample of nine such Milky Way-like 
haloes, which have Mh(2 = 0) = 2 x lO^^/i'^M©, where 
Mh is the total mass of the host halo, including the mass 
enclosed in dark matter substructures. In Fig. 1 we plot the 
evolution of their mass (upper- left panel), virial (lower- left 
panel) and scale radii (upper-right panel) over one Hubble 
time. As a comparison, the strong-solid line shows the av- 
eraged mass growth of a Milky Way-like halo obtained by 
van den Bosch (2002) using an extended Press-Schecter for- 
malism. In the lower-right panel, we also plot the variation 
of their concentration parameter (defined as the ratio be- 
tween the scale and virial radii) in time. It is interesting to 
note that both the mass and the virial radius of haloes grow 
almost linearly with time on average until the present epoch. 

For each realization in our sample, a different merger 
tree is defined. In Fig. 2 we show the number of substructures 
merging with the host halo per unit time (upper-left panel) 
and the cumulative number of merging events during one 
Hubble time (upper-right panel). Both panels show that the 
merging rate has slowed down since the early formation of 
the host halo. 

Each substructure has associated with it a virial mass. 
In the lower panels we plot the mass fraction /, which is 



t Benson (2004) finds some evidence for evolution of the velocity 
distribution to higher redshifts, and for dependence on the masses 
of the merging haloes. These dependencies are poorly quantified, 
and so we ignore them here. 
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Figure 1. Growth of the host halo in time. In each panel we plot 
the relevant property of each host halo from our sample of nine. 
Upper left-hand panel: The virial mass. Strong-full line shows the 
result of van den Bosch (2002) for a halo with Mh{z = 0) = 
2 X I0^^h~^ Mq Upper right-hand panel: The scale radius of the 
NFW density profile. Lower left-hand panel: The virial radius. 
Lower right-hand panel: The concentration parameter. 
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Figure 2. Upper-left panel: Evolution of the number of mergers 
per unit time. Upper-right panel: Cumulative number of merg- 
ers as a function of time. Lower-left panel: Evolution of the mass 

fraction of accreted substructures (see text) per unit time. Lower- 
right panel: Cumulative mass fraction as a function of time. 
Note that only substructures with an initial mass larger than 
10^ h~^ Mq were included in these calculations. 



defined as the mass enclosed in substructures, normalised to 
Mh{z = 0), which is accreted per unit time (lower- left panel) 
and the cumulative mass fraction (lower-right panel) as a 
function of time. This last quantity accounts for the mass 
gain of the host halo as a result of substructure mergers and, 
as we can see, ranges from 0.3-0.7 at z = 0. Interestingly, de- 
spite the decrease in the number of mergers, the cumulative 
mass fraction shows an approximately linear growth owing 
to the larger masses of the late-time merging substructures. 
For this plot, wc consider only those substructures with 
M > W^h~^MQ. We note that the mass fraction enclosed in 
subhaloes barely depends on the arbitrary minimum mass 
selected here since / is principally determined by a relatively 
low number of fairly massive substructures, in particular by 
those with M(i™) > 0.1Mh(z = 0) ~ 2 x 10"/i"^Mo (see 
Section 4.3.2), where tm is the merging time . 

In order to obtain a better understanding of the com- 
plex processes acting on dark matter substructures during 
their evolution we inject two control substructures in the ini- 
tial sets with orbital and internal properties at their merging 
time t = tm selected ad hoc. The scattering force induced 
by these bodies was settle to be zero for consistency. Con- 
trol substructures are assumed to follow NFW profiles. For 
these galEixies, we shall perform an illustrative study of the 
effects induced by substructure-substructure interactions. In 
Table 1 we show their orbital and structural parameters. 



3 THE SEMI-ANALYTICAL CODE 

In this section we present a description of the semi-analytical 

code used to model the evolution of the dark matter sub- 
structures described above as they orbit in the host poten- 
tial. 

3.1 Dynamics of dark matter substructures 

We treat dark matter substructures as point-mass particles 
that move within a host halo potential in order to compute 
their unperturbed orbits. Additionally, we include a dynam- 
ical friction force and account for forces between pairs of 
dark matter substructures by direct summation. 
Our model solves the equation of motion: 

f = fh + fsub + fdf, (1) 

where fh is force per unit mass (or specific force) induced 
by the host halo, fsub is that due to interacting dark matter 

substructures and fdf is the dynamical friction force exerted 
by the diffuse background of dark matter particles on the 
substructure. 

We assume that the host halo, as well as the dark mat- 
ter substructures, follow an NFW density profile (Navarro, 
Prenk & White 1997) 

pNFW = -7— I wT"^^ -, To' (■^) 
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rs 


^vir 


r(tm) 


e 






[kpc /i-i] 


[kpc /i-i] 


[kpc h-^ 




Control A 


108 


0.90 


6.24 


133.90 


0.00 


Control B 


108 


0.90 


6.24 


133.90 


0.50 



Table 1. Control substructures. Orbital and density parameters at t = tm- e denotes the orbital eccentricity, defined as e = (ra — 
'■p)/(^a + '■p), where ra and rp are the orbital apo and peri-centres, respectively. Control substructures are initially located at the 
apo-centre, r(tin). 



where po is a constant and rs is the halo scale length. The 
force, fh, acting on a particle at the position r can be ana- 
lytically derived 



fh = - 



GMh 



ln[l + r/rs] - r/(r + rs) 



(3) 



ln[l + rvir/vs] — rvir/ (rvir + rs) 

where Mh is the mass within the virial radius rvir after sub- 
tracting the mass bound in substructures, and r = r/r. Wo 
note that Mh, rs and rvir are time-dependent quantities (see 
Fig. 1). 

The presence of Ng^b substructures will perturb the or- 
bits that dark matter substructures would follow in a smooth 
potential. To estimate the perturbation we perform a simple 
force summation, so that the specific force acting on the i**" 
substructure is 



•GM, 



12)3/2 



(4) 



where is the softening parameter of the i"* substructure 
which accounts for its finite extent. Our approximation of 
point masses is valid providing encounters remain in the 
regime of \rj — r^l > r^^vir + rj.vir- Closer encounters, where 
two systems overlap, can be treated as inelastic collisions. 
In such large fraction of the orbital energy is trans- 

mitted to stars through resonances so that eq. (4) overes- 
timates the resulting velocity change. In order to compare 
our results to those of self-consistent N-body simulations we 
have chosen et = r^^h, ''i,h being the half-mass radius. This 
choice decreases the relative force for overlapping encounters 
and avoids a diverging force at very small relative distances, 
whereas if |rj — r^l > ri,vir + ^j.vir ~ 2ri,vir it recovers the 
Newtonian behaviour. 

Several studies of satellite orbital decay have shown 
that, in spherical systems, Chandrasekhax's formula for dy- 
namical friction (Chandrasekhar 1943) is sufficiently accu- 
rate if the Coulomb logarithm is treated as a free param- 
eter to fit to N-body data (e.g. van den Bosch et al. 1999, 
Colpi et al. 1999). Semi-analytic methods that include Chan- 
drasekhar's dynamical friction have been demonstrated to 
reproduce accurately the overall evolution of satellite galax- 
ies (e.g. Velazquez & White 1999, Taylor & Babul 2001, 
Benson et al. 2004) and, therefore, represent a useful tool 
for conducting extensive studies of a large parameter space. 
Chandrasekhar's formula can be written as 



fdf 



-47rGMp(< ■y)lnA-— ^ 



where InA is the Coulomb logarithm and 



p(< v) = pNFw(r) 



erf (A) —e 



(5) 



(6) 



where X = \w\/\/2ah and CTh is the one-dimensional ve- 
locity dispersion of the halo dark matter, defined as a = 

l/pNPw(r) pNFw(r')fh{r')dr' (e.g., eq. 2.12 of Hernquist 
1993). We fix In A = 2.1 from the results of Penarrubia, Just 
& Kroupa (2004). 

We use a Runge-Kutta algorithm to solve eq. (1) with a 
time-step At = 6.5 X W~*h~^ Gyr, short enough to properly 
resolve close encounters. Subsequently, we apply a leap-frog 
integration of the orbit at each time-step. 



3.2 Mass evolution 

Substructures experience mass loss whenever the external 
tidal field is sufficiently strong. Mass which becomes un- 
bound can subsequently escape from the satellite. Depend- 
ing on the rate of change of the external tidal field as seen 
along the satellite orbit, we cam distinguish between slow 
and rapid regimes: tidal limitation and tidal heating, re- 
spectively. 



3.2.1 Tidal limitation 

In the first regime, the unbound mass is determined by the 
tidal radius (King 1962), which is defined for a spherically 
symmetric satellite as the distance from the satellite centre 
at which the force due to the satellite's self-gravity and the 
external tidal force balance. If the satellite follows a circular 
orbit, the system can be considered as static in a rotating 
frame, and one can estimate the tidal radius as (King 1962) 



Rt 



a;2 - d2$h/dr2 



1/3 



(7) 



where ui is the angular velocity of the satellite and #h the 
potential of the host halo. 

This approximation is valid when: (i) The substructure 
mass is much smaller than that of the host halo and (ii) Rt 
is small compared to the distance to the centre of the host 
halo {Rt <^ r). Even under these conditions, the mass within 
Rt is not exactly equal to the bound mass, since there may 
be particles that have (small) positive energy and stay in 
transient orbits within the satellite (e.g. Binney & Tremaine 
1987). Furthermore, some particles at smaller radii, which 
nevertheless have orbits which take them to larger radii, can 
become unbound (Kampakoglou & Benson, in preparation). 

If the satellite follows a non-circular orbit one can still 
use eq. (7) to calculate the instantaneous tidal radius, with 
uj now being the instantaneous angular velocity. For these 
orbits, the mass loss is most rapid near peri-centre where 
the tidal force is greatest (e.g., Penarrubia, Kroupa & Boily 
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2002, Piatek & Pryor 1995). In our calculations, the effects 
of transient orbits are assumed to be negligible compared to 
the total amount of mass stripped out by the peri-galacticon 

passages. 

Whereas cq. (7) is quite accurate in accounting for the 
mass loss induced by tidal forces from the host halo, close en- 
counters between dark matter clumps may enhance the mass 
stripping process by heating (i.e. expanding) these systems. 
Such collisions occur on time-scales short compared to the 
orbital period and can be treated as tidal shocks. 



3.2.2 Tidal heating 

Satellites travelling through regions where the external po- 
tential changes rapidly suffer tidal shocks. This process can 
bo rnodollod as a perturbation, with a particular frequency 
spectrum, that adds energy to the satellite particles. 

Gnedin & Ostriker (1999) show that, as a result of the 
shock, the satcUite is contracted, with subsequent expan- 
sions and rc-contractions until it reaches a final state of 
equilibrium, in which the binding energy is reduced com- 
pared to the original value and the satellite has expanded. 
This non-equilibrium phase lasts, in their models, for around 
20 satellite dynamical times. 

The detailed study of tidal shocks is beyond the scope 
of this work. We choose to follow the method of Gnedin, 
Hernquist & Ostriker (1999), hereinafter GHO, in order to 
calculate the energy change of the substructure. 

In the impulse approximation, the mean change in the 
specific energies of particles in a spherical system, i, due to 
the encounter with an extended body, j, can be written as 



i(A.) 



Vob^ 



Rx{R) 



(8) 



where b, Vo are the relative distance and velocity between 
two substructures at the moment of closest approach, Mj is 
the mass of the system j and R is the particle's radius in the 
substructure. There are several corrections to this equation 
which must be taken into account: 

(i) Extended bodies: The function x(^) accounts for 

the fact that the energy change of particles belonging to the 
system i depends on the tidal field of the body j through its 
density profile. For the sake of brevity, we refer the reader 
to eq. (11)-(14) and eq. (17) of GHO for a description of the 
expression x(^)- 

(ii) Diffusion term: Gnedin & Ostriker (1999) show that 
the second-order energy gain (AE^) — ((vAv)'^) induces a 
shell expansion stronger than ( A_B) . This term can be imple- 
mented in our semi-analytic calculations as a multiplicative 
correction A to the first-order energy variation. Here we im- 
plement the results of Taylor & Babul (2001) who, using a 
similar mass-loss scheme, estimate A — 2.8 from fitting the 
mass evolution to N-body simulations. 

(iii) Adiabatic correction: In the impulse approxima- 
tion, one assumes that particles do not move during the en- 
counter. However, in reality particles in the substructure will 
move during the encounter, leading to a suppression of the 
heating, eq. (8) can be corrected for this fact by including 
an adiabatic correction term 



where x = tsh/toib{R), i.e., the ratio between the shock du- 
ration and the substructure's internal orbital period at a 
given radius. Specifically 

b_ 

R 



tah = 



(10) 



torb(-R) 



27r- 



iR)- 



and v(R) is the rotation curve of the substructure i. The ex- 
ponent 7 depends on the encounter duration relative to the 
half-mass dynamical time, tdyn of the system i, and varies 
from 7 = 2.5 for tsh < tdyn to 7 = 1.5 for tsh > tdyn (Gnedin 
& Ostriker 1999). The adiabatic correction tends to unity 
for X <^ 1 (typical of the outer parts of substructures) and 
to zero for a; ^ 1 (typical of the central parts). 

(iv) Time-VEirying correction: Eq. (8) uses the maxi- 
mum value of the tidal field (that is, the tidal field at the 
moment of closest approach between both substructures) in 
order to calculate the particles' energy change. However, in 
reality, the interaction occurs along the whole encounter. 
Gnedin & Ostriker (1999) show that the external tidal force 
can be approximated by a Gaussian centred at the moment 
of closest approach and with dispersion equal to the en- 
counter duration, obtaining a correction to eq. (8) from fits 
to N-body simulations that can be written as 



where y = tsh/Wn- 



(11) 



After combining eq. (8) with corrections (ii) through 
(iv), the resulting equation for the tidal heating induced by 
close encounters between substructures is 



AE = AEiXA{x)IcoTr{y). 



(12) 



To determine the changes experienced by the satellite 
after the energy injection, we assume that no shell cross- 
ing occurs. This approximation is quite accurate since the 
encounters between substructures lead to small changes of 
the internal energy due to their low masses and in this limit 
the change in radius of each shell, A_R, is a rapidly increas- 
ing function of radius (see below). Under this condition, a 
change of energy E —>■ E + AE results in an expansion of 
the system i 



AR. 



AER' 
GM,{Ry 



(13) 



^(a;) = (1-l-a;^ 



(9) 



so that the mass distribution after the tidal shock is 

M-{R') = M({R + AR) = Mi{R). As this equation suggests, 
some material will be expanded out of the tidal radius which 
will, therefore, enhance the degree of mass loss. Addition- 
ally, the mass loss induces a progressive shrink of the tidal 
radius since Rt oc M^^^. We note that the shell radii expand 
monotonically in time after each tidal shock. 

Our code computes the mass of substructure i at each 
time-step in the following way: 

1) Calculate Rt using eq. (7). 

2) If the distance between the substructure and any other 
substructure reaches a minimum during the time-step, the 
code computes AE and AR for each R in the substructure 
and subsequently computes the expansion of the substruc- 
ture. 



© 0000 RAS, MNRAS 000, 000-000 



6 J. Penarrubia & A. J. Benson 



3) Finally, the code calculates the new mass of the satellite 
M' = M(R < Rt). U M' > M then we force M' = M, so 
that the mass of the substructure decreases monotonically. 

As Taylor & Babul (2001) note, this technique suffers 
from some limitations: (i) The assumption of virial equilib- 
rium between the satellite shocks. (ii) Numerical simula- 
tions of tidal mass stripping have shown that the process is 
highly anisotropic. Our mass loss scheme docs not account 
for any possible net angular momentum of the escaping par- 
ticles, (iii) Overlapping encounters cannot be accounted for 
by this simple mass-loss scheme. In such cases, depending 
on the mass ratio of the clumps and their relative velocity, 
the encounter may result in the destruction of one or both 
substructures. Even taking account of these limitations, this 
scheme has been proven to lead to accurate estimates of the 
mass evolution (Taylor & Babul 2001, Penarrubia 2003). 



4 RESULTS 

The goal of this paper is to develop a better understanding 
of the dynamical processes that may influence the present 
distribution of dark matter substructures. The main advan- 
tage of the semi-analytic code outlined above is that we can 
study separately the effects induced by each process. 

In Table 2 we describe the calculations carried out for 
this work. Simulations C1,C2,C3 and C4 combine the imple- 
mentation of dynamical friction and gravitational scattering 
in the semi-analytic code. The calculation C5 switches off 
the mass evolution scheme by imposing a constant mass for 
all substructures during their evolution. Lastly, simulation 
C6 repeats the calculation of CI, but with a minimum mass 
at merging for substructures which is ten times larger than 
for the other simulations. 

For each of the calculations, we evolve the substructure 
samples described in Section 2. 

4.1 Interactions between dark matter 
substructures 

The internal evolution of dark-matter substructures might 
be affected by repeated encounters with other substructures. 
In contrast with the effects of two-body collisions with par- 
ticles of the smooth dark matter component (which can 
be treated analytically as dynamical friction) , substructure- 
substructure encounters may lead to a strong tidal heating 
(i.e. to an overall expansion of the system) which would ac- 
celerate mciss stripping due to tidal forces. That physical 
process, however, will induce significant effects only when 
the encounter is sufficiently close. 

In order to illustrate the properties of two-body encoun- 
ters between substructures within a given host halo, we use 

i The estimation of the mass shell expansion apphed in this work 
assumes that an extended system virializes after a tidal shoek 
event. Gnedin & Ostrikcr (1999) show that the viralization pro- 
cess lasts approximately 20 dynamical times of the system. Thus, 
our method may not consistently reproduce the subhalo heating 
if the system suffers strong encounters within short time intervals 
(fortunately, we shall show in Section 4.1 that such intera<;tions 
are fairly rare). 



two substructures (control A and B, see Table 1) injected in 
each of our host halo samples with kinematical and internal 
properties selected ad hoc. The full dynamical scheme (CI) 
was selected. For each of those substructures, we calculate 
at each time step the relative distance to all other substruc- 
tures. If a minimum occurs, we consider that an encounter 
takes place and record the orbital parameters. 

In Fig. 3 we plot the distribution of encounter param- 
eters for all encounters that our control substructures ex- 
perience during the last 5.5 h^^Gyv in all host haloes. In 
the upper-left panel we show the distribution of minimum 
encounter distances for our control substructure A (circular 
orbit) and B (eccentric orbit, e — 0.5). Both substructures 
sfiow curves with a marked maximum at 6 ~ 85h^^ kpc 
(comparable to the size of the haloes). For both orbits, the 
number of close encounters 6 « is fairly small, as can be 
seen from the upper-right panel, so that during the last 7 
Gyr the average number of overlapping encounters (i.e. those 
with minimum encounter distances of feJi 2rvir — lOh^^ kpc, 
where rvir is the substructure's virial radius) was less than 
10, independently of their initial eccentricity. This result 
favours our semi-analytic calculation of substructure orbits 
since relative forces along non-overlapping collisions can be 
accurately reproduced by interactions of point-masses. It is 
also interesting to note that the substructure following an 
eccentric orbit has suffered around 33% more encounters 
than that following a circular one. The distribution of rela- 
tive velocities at the moment of closest approach also shows 
unique maxima at Vo ^ 320 km/s (control A) and Vo ^ 340 
km/s (control B). Lastly, in the lower-right panel we plot 
the evolution of the encounter rate dN/dt for both control 
substructures. Control A shows a linear growth of dN/dt, 
reflecting a progressive increase in the number of substruc- 
tures within the host galaxy (see Fig. 2). The oscillations 
in the encounter rate of Control B are due to its eccentric 
orbit and show that the number of encounters is enhanced 
at the peri-centre of the orbit. As for control A, the mean 
rate within one orbital period grows linearly in time. 

Although the results obtained in this section are purely 
illustrative, we may foresee that overlapping collisions are 
rare, representing less than 2% of the total number of colli- 
sions for our control substructures. More distant encounters 
will induce negligible tidal heating since, from eq. (8), we 
have that the heating energy scales as AEi oc 1/6*. 



4.2 Evolution of dark matter substructures 

In this section we analyse the effects that substructure- 
substructure interactions induce on the final distribution of 
substructures within the host halo. 

We note that, in order to compare our results to those 
of recent N-body calculations, our statistical study was car- 
ried out in a way similar to that used when analysing cos- 
mological N-body simulations: (i) each halo realization was 
evolved separately; (ii) the resulting substructure samples 
were compiled into a single sample; (iii) we only consider 
substructures within the virial radius of each of the host 
haloes and with final mass larger than the resolution mass 
(which we introduce artificially into our calculation for the 
purposes of this comparison) . 
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Name mass resol. (m) mass loss dyn. friction scattering 

[/i-IMq] 



CI lO'^ ("^2) yes yes yes 

C2 lO'' (rrt2) yes no yes 

C3 10^ ("^2) yes yes no 

C4 10^ ("12) yes no no 

C5 10^ ('"■2) no yes yes 

C6 10* (mi) yes yes yes 



Table 2. Calculations carried out. Wc study different mass resolutions by evolving only those substructures that at the time of merging 
have M > m. We also analyse the effects that dynamical friction and gravitational scattering induce on the evolution of substructures. 
We note that in those calculations without substructure scattering, the tidal heating scheme was also switched off for consistency. 




V„ (km s"') t (Gyr h"') 



Figure 3. Distributions of minimum approach distances (upper- 
left panel) and relative velocities at the moment of minimum ap- 
proach (bottom-left panel) for the control galaxies A (solid lines) 
and B (dotted lines) of Table 1. In the upper-right panel we also 
plot the cumulative number of encounters as a function of mini- 
mum distance. In the bottom-right panel we show the encounter 
rate as a function of time. 



4-2.1 Structural evolution 

Following the above prescriptions, we determine, for the cal- 
culation CI, the average substructure expansion (heating) 
as a consequence of substructure scattering. In Fig. 4 we 
show the fraction of substructures as a function of the ex- 
pansion ratio of the scale-length (full line) and the virial 
radius (dotted line). As can be seen, the fraction of sub- 
structures that were inflated by more than 10% of their ini- 
tial size (logio[r(t/)/r(tm) — 1] > —1) due to tidal heating 
corresponds to approximately 2% of the total number of sub- 



Figure 4. Fraction of substructures as a function of their expan- 
sion ratio. The full and the dotted lines show the expansion ratio 
of the scale-length (r^) and the virial radius (rvir), respectively. 

structures. The remaining substructures suffer an expansion 
that can be considered negligible for mass loss computations. 

In light of this result, we expect tidal heating to induce 
negligible effects on the final substructure mass function. 

4-2.2 Orbital evolution 

In this work we have assumed that the interaction of sub- 
structures with the rest of the dark matter present in the 
host halo can be separately treated as: (i) dynamical fric- 
tion (smooth dark matter component) and (ii) substructure- 
substructure scattering (clumpy dark matter component). 
We have analysed the effects that each component induces 
on the orbits of substructures by switching on and off these 
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terms in eq. (1). Due to the growth in mass and size of the 
host halo, adiabatic invariants (e.g see Section 3.6 of Bin- 
ney & Tremaine 1986 for details) are adequate quantities to 
reflect how those processes alter the orbits of substructures 
since thoy remain constant during adiabatic (slow) changes 
of the host halo potential. Thus, in our calculations, their 
values can only be altered either by dynamical friction or by 
substructure-substructure interactions. 

The total angular momentum of systems orbiting in 
isotropic potentials is an adiabatic invariant. If only dy- 
namical friction is implemented, the angular momentum de- 
creases as a consequence of the drag force, whereas if only 
scattering is implemented substructures can gain or lose an- 
gular momentum depending on the characteristics of the 
two-body encounters occurring during their evolution. In 
Fig. 5 we plot the fraction of substructures as a function of 
their angular momentum variation. The dotted line, calcula- 
tion C2 (i.e. switching on dynamical friction) shows that the 
effect of the smooth dark matter component is practically 
negligible, since fewer than 1% of the surviving substruc- 
tures have suffered a measurable angular momentum loss. 
This result is fully in agreement with that of Zhao (2004,) 
who, on the basis of semi-analytical calculations of satel- 
lite orbits in a growing host potential, conclude that the 
effects of dynamical friction are suppressed by mass loss. 
On the contrary, the effects of substructure-substructure in- 
teractions (calculation C3; dashed line) on individual orbits 
of dark matter clumps are considerably stronger. The distri- 
bution of the angular momentum variation follows approx- 
imately a Gaussian centred at L{tf)/L{tm) ~ 1 with dis- 
persion of around 10%. This figure shows that substructure 
scattering has altered the orbits of the present substructures 
much more effectively than dynamical friction. Implement- 
ing both dynamical friction and scattering (calculation CI) 
produces basically no further change in the angular momen- 
tum distribution. 

On the basis of these results, the following question 
arises: Why are dynamical friction and gravitational scat- 
tering inefficient at altering substructure orbits? To answer 
that question we compare in Fig. 6 the averaged angular mo- 
mentum change {L(t f) / L{tm)) and the deviation from the 
average a (i.e. the variance of the distribution) of substruc- 
tures evolved in calculations where our mass loss scheme has 
been switched on (full line; calculation CI) and off (dotted 
line; calculation C5). We plot those quantities as a function 
of the substructure initial mass (i.e. at the merging time) 
and the present orbital radius. 

Looking first at the dotted lines, we observe that if 
we allow substructures to evolve until 2 = (i.e. they are 
not disrupted by tidal fields), the averaged angular momen- 
tum change is considerable, independent of initial mass and 
present radius. Moreover, the angular momentum variation 
as a function of mass (lower-left panel) resembles that typ- 
ically found in stellar systems with a low number of par- 
ticles and a given mass spectrum where, owing to two- 
body relaxation processes, the angular momentum of most 
massive bodies decreases on average (see, e.g.. Section 8.4 
of Binney & Tremaine 1986). This panel shows that sub- 
haloes with M(tm) > 0.01Mh(z = 0) lose angular momen- 
tum {{L\tf\/ L\tm\) < 1) on average whereas light subhaloes 
M{tJ) < 0.01Mh(^; = 0) gain it {{L[tf\/L[tm\) > 1). The 
lower-right panel also shows that the deviation from the av- 



-1 



CI 
C2 
C3 



-0.2 



0.2 



log,„[L(t,)/L(t„)] 



Figure 5. Fraction of substructures as a function of the angu- 
laj: momentum variation. We show the distribution resulting from 
calculations where (i) the force term accounts for dynamical fric- 
tion and substructure scattering (full line), (ii) the dynamical fric- 
tion term is switched off (dashed line) and (iii) the substructure- 
substructure scattering term is switched off (dotted line). Sub- 
structures included in this sample have M > 10^ Mq and are 
located at r < rvir at z = 0. 



erage (a) is also larger for substructures that have lost a 
large fraction of their initial angular momentum. 
Besides the dependence of the angular momentum distribu- 
tion on the initial mass of substructures we also find a radial 
dependence. In the upper-left panel we show the averaged 
angular momentum variation as a function of radius. The 
dotted line shows that substructures located at r < 0.12rvir 
lose on average approximately 60% of their initial angular 
momentum, whereas the angular momentum of substruc- 
tures at r > O.SSrvir has increased by 10-20%. From these 
results, it appears that changes in angular momentum and 
orbital radius are related, with massive substructures decay- 
ing to the inner regions of the host halo, while light substruc- 
tures are ejected to the outer regions as a result of two-body 
interactions. 

The present angular momentum distribution is strongly 
modified if we allow substructures to lose mass (full line). 
In that case, {L{t f) / L{tm)) — 1, independently of mass 
and radius. Only the most massive substructures (M[im] > 
0.05Mh[2; = 0] ~ 10^^ft~^MQ) show a decrease of angular 
momentum (up to 28%). Substructures at r < 0.2rvir also 
present some degree of dynamical evolution, with a distribu- 
tion centred at {L{tf)/L{tm)) — 1 but with a considerably 
larger deviation from the mean value cr^ ~ 0.07 (in contrast 
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Figure 6. Averaged change of angular momentum as a function 
of radius (upper- left panel) and initial mass (lower-left panel) . We 
also plot the variance of the distributions as a function of radius 
(upper-right panel, denoted as ar) and initial mass (lower-right 
panel, denoted as um)- We compare calculations where mass loss 
is switched on (full line; CI) and off (dotted line; C5). 



with those substructures at r > 0.2rvir, which show ar ^ 
0.01-0.02). 

These results appear to indicate that the meiss loss pro- 
cess prevents significant alterations of the initial adiabatic 
invariants, thus minimising dynamical evolution effects such 
as orbital decay and gravitational scattering. Only massive 
subhaloes (M[tin] > 0.05Mh[2: = 0]) show a measurable or- 
bital evolution, but these represent less than 2% of the total 
number of substructures in our sample. 



4.3 Present distribution of substructures 

4-3.1 Spatial distribution 

In this Section we analyse how dynamical friction and sub- 
structure scattering affect the final radial distribution of 
substructures. In order to compare our results with those 
from cosmological calculations at z = 0, we mimic the cri- 
teria used to identifying substructures: our subhalo sample 
includes only those substructures with masses above a cer- 
tain value (i.e mass resolution). In N-body calculations, that 
threshold is determined by the number of particles and the 
mass of the host system since only those systems with a 
minimum number of particles (equivalently, mass) can be 
distinguished as bound substructures at any time of the sim- 
ulation. Additionally, we investigate the effects of the mass 



resolution by performing calculations with two different min- 
imum masses (similar to those of Gao et al 2004). 

In the upper panel of Fig. 7 we plot the present radial 
distribution of substructures for calculations where the dif- 
ferent force terms of oq. (1) where switched on and off (see 
Table 2 for details). Additionally, the dot-dashed line shows 
the resulting distribution of sustructures if the mass loss 
scheme is not implemented and the heavy solid line shows 
the mass distribution of the host halo. This panel shows 
that: 

(i) The cumulative profile of substructures is on all scales 
less concentrated than that of the smooth dark matter com- 
ponent. 

(ii) The largest difference in the final distribution is induced 
by mass loss. Switching on our mass loss treatment drasti- 
cally reduces the number of substructures on all scales. 

(iii) The dynamical friction and scattering terms of eq. (1) 
ordy introduce differences of around 5% in the number of 
subhaloes at r < O-SSrvir. 

(iv) Our semi-analytical treatment allows us to disentangle 
the contribution of dynamical friction and that of scatter- 
ing. A comparison between the calculations 02 (scattering)- 
03 (dynamical friction) shows that gravitational scattering 
slightly increases the number of subhaloes at r < 0.35rvir, 
whereas if comparing 03-04 (no dynamical friction, no scat- 
tering) we observe that the effects of dynamical friction can 
be neglected at all radii. 

In the lower panel of Fig. 7 we have analysed the ef- 
fects of the mass resolution on the spatial distribution of 
substructures. We have caxried out calculations where the 
minimum initial mass of merging substructures was above 
two different thresholds mi = lO^h-^MQ (calculation 06) 
and m2 = 10^ Mq (calculation 01); in other words, we 
vary the substructure mass resolution, denoted here as mi, 
m2, where m2 = mi/10. For determining the present dis- 
tribution we also choose those substructures with masses 
above mi (dotted line) and m2 (full line) located within 
the virial radius. The resulting cumulative profiles are com- 
pared against the fitting expression obtained by Gao et al. 
(2004) from large cosmological N-body simulations, finding 
that the semi-analytic curves presented here agree extremely 
well with the N-body results. Only for r < 0.2rvir does the 
abundance obtained semi-analytically become smaller than 
that predicted by the formula of those authors. However, 
this problem is likely due to their selection of the fitting ex- 
pression, since a similar mismatch also appears in their fit 
(see their Fig. 11). As Gao et al. (2004) claim, the mass reso- 
lution scarcely affects the distribution of substructures, only 
in the innermost region of the host halo is the abundance 
of substructures slightly suppressed when reducing the mass 
resolution. 



4.3.2 Mass function 

We repeat the analysis of Section 4.3 for the substructure 
mass function. In Fig. 8 we show in the upper-left panel 
the curves resulting from calculations where the interaction 
forces (i.e. dynamical friction and substructure-substructure 
scattering) were switched off (dotted line, calculation 04) 
and on (full line, calculation 01). As we can see, the dif- 
ferences are practically negligible on all mass scales. The 
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Figure 7. Cumulative radial distribution of substructures at 
z = 0. Upper panel presents the distribution resulting from using 
different force terms and that obtained if the mass loss treatment 
is switched off (see Table 2 for details). We also show for compar- 
ison the cumulative profile of the host halo. We note that those 
substructures included in the distribution have M > lO"^ h^^ Mq 
at 2 = and are located within the virial radius. In the lower-left 
panel we compare the cumulative profile that results from simula- 
tions where the minimum mass of substructures at their merging 
time was above lO^ft-^M© (dotted line) and lO^/i-iM© (full 
line). The strong full line shows the results of Gao et al. (2004) 
from large N-body simulations. 

dashed linos shows the resulting mass function if keeping 
constant the mass of substructures during their evolution 
(calculation C5). In that case, substructures as massive as 
M = 0.25Mh(z = 0) = 5 X lO" /i^^M© can be found in 
the final sample, whereas the total number of substructures 
is around a factor 3 larger than in calculations where the 
mass loss scheme is implemented. Moreover, it is remark- 
able that both curves present a similar slope, which is due 
to the fact that, in average, all substructures lose a similar 
mass fraction (see Section 4.4). 

In the upper-right panel, we plot the averaged ratio be- 
tween the mass enclosed in substructures and the host halo 
mass (which we denote as (/)) as a function of the substruc- 
ture mass. This panel shows that the present population of 
substructures represents only 2.8% of the mass of a Milky 
Way-like halo. This result is in very good agreement with 
that of van den Bosch, Tormen & Giocoli (2004), who find 
(/> = 2.5% for haloes with Mi,{z = 0) = 2 x W^^h'^MQ. 
Wo find that switching off mass loss loads to a mass frac- 
tion of {/) ~ 25%, so that the mass lost by substructures 
contributes in a signiflcamt amount to the host halo mass. 
We note that the averaged mass fraction found here is con- 



siderably smaller than that plotted in the lower-right panel 
of Fig. 2 (which results in (/) ~ 0.44), since for this calcu- 
lation we excluded all substructures at r > rvir. Like van 
den Bosch, Tormen & Giocoli (2004), we find that most of 
the present day mass fraction originates from substructures 
with M/Mh > 10"'^ if implementing mass loss (note that the 
curves are remarkably flat for lower mass values). However, 
the largest contribution is shifted to M/Mh > 10~^ if no 
mass loss occurs, so that the mass enclosed in substructures 
with M/Mh < 10^^ represents a tiny fraction of the entire 
mass despite their larger abundance. 

In the lower panels, the above analysis is repeated for 
calculations with a mass resolution ten times smaller. Ad- 
ditionally, we reproduce the fornmla suggested by van den 
Bosch et al. (2004) (strong full line) to fit the mass functions 
obtained from cosmological N-body calculations. As in the 
above panels, we obtain practically no difference in the final 
mass function. 

Lastly, we must remark the good agreement between the 
N-body mass function and that obtained semi-analytically 
here§ . 

4.4 Mass evolution of substructures 

In this Section we analyse the effects that dynamical friction 
and substructure scattering induce on the mass evolution of 
dark matter substructures. Additionally, we also determine 
the dependence of the results on the mass resolution. 

In Fig. 9 we plot the fraction of substructures as a func- 
tion of their mass decrease (left column) and the averaged 
mass loss as a function of the substructure present mass 
(right column). In the upper row we analyse the effects of 
dynamical evolution, whereas in the lower row we compare 
the results from calculations with different mass resolution. 
In the upper-left panel we can see that dynamical friction 
and gravitational scattering induce slight changes on the 
mass loss process. The distribution peaks at M{tf) ~ 0.5- 
0.8M(tm). Strong mass loss events, where m2 < M{tf ) < 
10^'^ M{tin) appear fairly rare. It is also interesting to note 
that a large number of substructures have lost a negligible 
fraction of their initial mass, either because they have been 
accreted relatively late or because their orbits keep them 
in the outer region of the halo, where tidal fields are weak. 
In the lower-left panel we can also see that mass resolution 
scarcely affects the mass fraction lost by substructures dur- 
ing their evolution. 

In order to infer a possible mass dependence in the mass loss 
process, we plot the averaged M{tf)/M{tin) as a function 
of the present mass (right column). Interestingly, we find 
that substructures lose, on average, around 60% of their ini- 
tial mass for M{tf) < 10"^Mh(t/), and 40% for those with 
M{tf) > W~^Mh{tf) ,although that difference might be 

§ The over-prediction in the number of substructures with M > 
O.lMij is a result of the universality of the fit (in the sense that 
they do not include the host halo mass dependence) carried out 
by Gao et al. (2004). Numerical and semi-analytical calculations 
of mass functions (e.g van dcii Bosch et al. 2004 and references 
therein) show that, for Mh(2 = 0) = 2 x lO^'^h'^MQ, a high- 
mass cut-off appears at Mcut — lO^^?i~^M0, very similar to that 
found here (Mcut = 1-2 x lO" /i-^Mq). 
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Figure 8. Left column: Substructure mass function at 2 = 0. 
Riglit column: Cumulative mass enclosed in substructures nor- 
malised to the host halo mass (for brevity, mass fraction). Upper- 
left panel shows the curves resulting from the implementation of 
dynamical friction and substructure scattering in eq. (1) (full line) 
and those where those forces were ignored (dotted line). For com- 
parison, the dashed line accounts for the initial mass function 
(i.e. switching off mass loss). In the lower-left panel we plot the 
mass function for two different mass resolutions mi = 10*/i~^Mq 
(dotted line) and 1712 = lO'' h~^MQ (full line). 

caused by the relatively small number of massive substruc- 
tures at « = (note that they represent less than 1% of 
the total number of substructures in our sample, see Fig. 8). 
This result is mostly independent of dynamical evolution 
processes and mass resolution. 



5 CONCLUSIONS 

In this paper, we have carried out semi-analytic calculations 
of the evolution of dark matter substructures in a Milky 
Way-like halo. Wc have generated nine samples of substruc- 
tures with spherical symmetry by following the technique of 
Cole et al. (2000) , which constructs the merger history of a 
halo back to high redshift. Subsequently, our semi-analytical 
code calculates the orbit of each substructure from the time 
of merging until the present epoch by solving its equation 
of motion. Additionally, this code implements a mass loss 
scheme that has been tested in previous papers and in- 
cludes an analytical treatment of the tidal heating induced 
by substructure-substructure collisions. 

We first use our code to carry out an illustrative exam- 
ple that quantifies the frequency and severity of two-body 
scatterings. We inject into each of our samples two control 



Figure 9. Left column: Fraction of substructures as a func- 
tion of the ratio of present mass to mass at time of accretion 
{M{tf)/M(tjn))- Right column: Averaged mass loss as a function 
of the present mass M(tf). Upper raw; Effects of dynamical evolu- 
tion on mass loss. We compare the results from calculations with 
and without dynamical friction and gravitational scattering (full 
and dotted lines, respectively). Lower raw: Effects of mass resolu- 
tion. We compare the results from calculations with mass resolu- 
tion 1712 = Wh-^MQ (full lines) and mi = lO^/i-iM© (dotted 
lines). For this plot, only substructures with masses above the 
mass resolution (either mi or m2, depending on the calculation) 
were included in the distribution. 

substructures that initially follow a circular and an eccentric 
orbit. The encounter parameters are recorded and statisti- 
cally analysed. We find that overlapping encounters (those 
where the minimum distance is less than 2rvir) are fairly 
rare, representing less than 2% of the total number of im- 
pacts. As a consequence, the tidal heating experienced by 
substructures owing to substructure-substructure encoun- 
ters induces a negligible expansion of their mass profiles. 

The main purpose of our semi-analytical investigation 
is to disentangle the effects of those processes acting on the 
evolution of substructures, which would be extremely expen- 
sive in terms of computational time if using large cosmolog- 
ical N-body calculations. 

After reproducing the most recent N-body results, we have 
separately analysed the forces that mainly drive the orbits 
of substructures, which are: (i) the force induced by the host 
halo potential, (ii) dynamical friction and (iii) substructure- 
substructure scattering. Additionally, we have also varied 
the mass resolution of our calculations by selecting only 
those substructures with a mass at their merging time above 
two different minimum masses lO^ft-^Mo and lO^h'^MQ. 
We find that: 
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(i) Neither dynamical friction nor substructure-substructure 
scattering induce strong variations on the present radial 
profile of substructures. Gravitational scattering slightly in- 
creases the number of substructures in the inner regions of 
the host halo, r < 0.35rvir, whereas the effects of dynamical 
friction can be neglected at all radii. The maximum increase 
in the cumulative number of substructures is less than 5% 
if compared with the distribution where both terms were 
switched off. 

(ii) The eflfects of those forces on the present mass function 
are negligible. 

We have studied the orbital evolution of substructures 

by measuring their angular momentum each time-step. In 
contrast to the energy, the angular momentum is an adia- 
batic invariant that remains constant in varying potentials 
in the absence of external forces. We find that the most mas- 
sive subhaloes suffer some degree of angular momentum loss, 
{L{z = 0)/L(tm)) ^ 0.7 for M{z = 0) ~ 0.1Mh(^ = 0). For 
less massive substructures, the averaged angular momentum 
remains practically unaltered since the time of their merg- 
ing. 

We have also determined how the mass evolution of sub- 
structures affects their present distribution. Wo find that 
mass loss effects are much more important than those in- 
duced by dynamical friction and scattering. Mass loss pro- 
cesses decrease the maximum mass expected for substruc- 
tures at the present day from Afcut ~ 0.25Mh down to 
Afcut = 0.06Mh — values in fairly good agroomont with those 
reported by Gao et al. (2004). We find that the mass con- 
tribution of disrupted substructures to the present day halo 
mass is around 22%, whereas only 0.025Mh is still enclosed 
in form of bound substructures. The present mass fraction 
in subhaloes mostly originates from a relatively low number 
of massive substructures (M ~ 0.1Mh[2 = 0]). As a conse- 
quence, the mass fraction is barely dependent on the mass 
resolution. These results agree fairly well with those of van 
den Bosch et al. (2004). Lastly we find that, on average, sub- 
structures lose around 40-60% of their mass by the present 
day, independent of their initial mass. This explains the low 
sensitivity of the mass function slope to mass loss processes. 

We conclude, therefore, that the present distribution 
of substructures is mainly determined by the merger tree 
and the mass evolution of these systems, whereas dynamical 
friction and gravitational scattering play a negligible role. 

These results are also important for understanding the 
present mismatch between CDM predictions and observa- 
tions of satellite galaxies: 

(i) The missing satellite problem. The apparent over- 
estimation in the number of satellite galaxies (o.g Klypin, 
Kravtsov & Valcnzucla 1999) cannot be explained by evo- 
lutionary processes within the host galaxy. The suggestion 
that dynamical friction or scattering may strongly enhance 
the disruption of substructures by driving them to the inner 
regions of the galaxy can be discarded. 

(i) The Holmberg effect. The anisotropic distribution of 
satellite galaxies (Holmberg 1969, Zaritsky & Gonzalez 1999, 
Sales & Lambas 2004), where most appear to follow nearly 
polar orbits out to 500 kpc, cannot be explained by dynam- 
ical processes acting within the host halo. 

Another interesting consequence of our results is that 
satellite galaxies beyond CSSrvir — WOh~^ kpc with masses 
M{z = 0) < O.OSMh ~ 10"/i"^M© have adiabatic invari- 



ants that have not been altered by dynamical processes act- 
ing within the host galaxy. Thus, these objects retain their 
orbital properties which are determined a priori by the in- 
tergalactic environment. Although we must keep in mind 
that observational surveys arc usually biased towards mas- 
sive satellite galEixies, which may present some degree of 
dynamical evolution (in particular we find that subhaloes 
with M{z = 0) > O.OSMh show a decrease of their initial 
angular momentum of around 25-30% on average) , observer 
tions of satellite orbits should provide a useful tool to infer 
the intergalactic dark matter distribution. 
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